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A heavy top with a fixed point and a rigid body in an ideal fluid are important 
t — ■ • examples of Hamiltonian systems on a dual to the semidirect product Lie algebra 

e(n) = so(n) k R™. We give a Lagrangian derivation of the corresponding equations of 
motion, and introduce discrete time analogs of two integrable cases of these systems: 
the Lagrange top and the Clebsch case, respectively. The construction of discretiza- 
tions is based on the discrete time Lagrangian mechanics on Lie groups, accompanied 
by the discrete time Lagrangian reduction. The resulting explicit maps on e*(n) are 
Poisson with respect to the Lie-Poisson bracket, and are also completely integrable. 
Lax representations of these maps are also found. 



Abstract 



(N 

> ! 1 Introduction 

<N 

The rigid body dynamics are rich with problems interesting from the mathematical point 
' of view, in particular, with integrable problems. Certainly, the most famous ones are 

the three integrable cases of the rotation of a heavy rigid body around a fixed point in 
a homogeneous gravity field, named after Euler, Lagrange, and Kovalevskaya. Another 
problem with a number of integrable cases is the Kirchhoff 's one, dealing with the motion 
of a rigid body in an ideal fluid. Here integrable cases carry the names of Kirchhoff, 
Clebsch, Steklov and Lyapunov. All these classical cases were discovered in the 18th and 
j> \ the 19th century. However, the list of intergrable problems of the rigid body dynamics is 

by far not exhausted by these ones. For example, the integrablity of equations describing 
the rotation of a rigid body around its fixed center of mass in an arbitrary quadratic 
potential is a much more recent observation due to Reyman p9l and Bogoyavlensky |T 



(However, some particular case of this result was given already by Brun [12|; the equations 
of motion in this case are identical with those describing the integrable case of the motion 
of a rigid body in an ideal fluid, due to Clebsch jnj). 

It was realized already in the 19th century that some of these problems admit multi- 



dimensional generalizations. According to bibliographical remarks in [ 10 1 and in |TH], the 



multi-dimensional generalization of the Euler top is due to Frahm [16| and Schottky [32]. 



A modern period of the interest in this system started with the general theory of geodesic 
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flows on Lie groups, due to Arnold ||, as well as with the paper by Manakov 2C], who 
found the n-dimensional Euler top as a reduction of the so called n wave equations, along 
with a spectral parameter dependent Lax representation. This was followed by a cycle of 
papers by Mishchenko and Fomenko and their school, generalizing the construction of a 
rigid body for a wide class of Lie algebras, see a review of this work in [36]. 

As for the multi-dimensional analogs of the Lagrange top, there exist two different 
versions. The first one, living on so(ra) >< W 1 , was proposed in [29[| , see also [31]. 
The second one, living on so(n) x so(n), was introduced in [27]. The multi-dimensional 
generalization of the Kovalevskaya top is due to |3Q| , see also [31, |j. For the Clebsch 
problem, the multi-dimensional generalization was found in |2q1 , for the Steklov-Lyapunov 
case a multi-dimensional generalization, living on so(n) ix so(n) and not on e(n) = so(n) x 
W 1 , was found in 0]. 

In the present paper, we find integrable discretizations for the multi-dimensional La- 
grange top and for a certain flow of the integrable hierarchy of the multi-dimensional 
Clebsch case. Our construction is based on the discrete time Lagrangian mechanics, as 
pioneered by Moser and Veselov [25], and further developed in ||, 34]. In order to give 
the reader a flavour of our integrable discretizations, we present here the corresponding 
formulas. In these formulas (M, P) £ e*(n), i.e. M € so(ra), P S W 1 . The equations 
of motion of the totally symmetric Lagrange top in the moving frame (which are also 
equivalent to the equations of motion of an arbitrary Lagrange top in the rest frame) and 
their integrable discretization are (^4 € W 1 is a constant vector): 




M k+l 
Pk+i 



--M k + eAAP k , 

I - (e/2)M fc+1 
I+(e/2)M k+1 



P k . 



Similarly, the equations of the Clebsch case of the rigid body motion in an ideal fluid and 
their integrable discretization read (B = diag (b\, . . . , b n ) is a constant matrix): 



M = PA(BP), 
P = -MP, 



Mi 



k+l 



k+1 



--M k + e 
I 



Pk A (BP k 



i + ( e yA)(p k ,BP k ) 

(e/2)Af fc+ i 



I+(e/2)M k+1 



Pk. 



We recall the general theory of the Lagrangian reduction in the continuous time and in 
the discrete time contexts, respectively, in Section 2 and 3. Further, we give in Section 4 a 
Lagrangian derivation of the equations of motion of a multi-dimensional heavy top, and in 
Section 5 we discuss the Lagrangian theory of the multi-dimensional Lagrange top, both 
in the moving frame and in the rest frame. Section 6 is devoted to the construction of an 
integrable Lagrangian discretization of the multi-dimensional Lagrange top, again in both 
formulations (in the moving frame and in the rest frame). Finally, in Section 7 we discuss 
the Lagrangian formulation of the problem of the rigid body motion in an ideal fluid, its 
integrable case due to Clebsch, and construct an integrable Lagrangian discretization of 
this problem. Conclusions are contained in Section 8. 
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2 Lagrangian mechanics and Lagrangian reduction on TG 

Recall that a continuous time Lagrangian system is denned by a smooth function L(g, g) : 
TG h M on the tangent bundle of a smooth manifold G. The function L is called the 
Lagrange function. We will be dealing here only with the case when G carries an additional 
structure of a Lie group, with the Lie algebra q. For an arbitrary function g(t) : [to, t\] i— > G 
one can consider the action functional 

h 

L(g{t),g(t))dt. 

to 

A standard argument shows that the functions g(t) yielding extrema of this functional (in 
the class of variations preserving g(to) and g(ti)), satisfy with necessity the Euler-Lagrange 
equations. In local coordinates {g 1 } on G they read: 

d /dL\ _ <9L 

dt \WJ ~w 

The action functional S is independent on the choice of local coordinates, and thus the 
Euler-Lagrange equations are actually coordinate independent as well. For a coordinate- 
free description in the language of differential geometry, see Q |22| . 
Introducing the quantities 

n = VgL e t;g, 

one defines the Legendre transformation: 

(g,g)eTG -> (g,U) eT*G. 

If it is invertible, i.e. if g can be expressed through (g, LT), then the Legendre transformation 
of the Euler-Lagrange equations (2.1) yield a Hamiltonian system on T*G with respect to 
the standard symplectic structure on T*G and with the Hamilton function 

H(g,U) = (U,g)-L(g,g), (2.2) 

(where, of course, g has to be expressed through (g, II)). 

When working with the tangent bundle of a Lie group, it is convenient to trivialize it, 
translating all vectors to the group unit by left or right multiplication. We consider first 
the left trivialization: 

(g,g)eTG » (j,!l)eGx e , (2.3) 

where 

Q = L g -i*g <3> g = L gif VL. 

Denote the Lagrange function pushed through the above map, by L^(g, Vl) : G x q h- > R, 
so that 

L^(g,n)=L(g,g), 

The trivialization (^^) of the tangent bundle TG induces the following trivialization 
of the cotangent bundle T*G: 



(g,U)GT*G ^ (g, M) £ G x g*, 
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where 



M = L*U o- n = L*_ 1 M. 



We now consider the Lagrangian reduction procedure, in the case when the Lagrange 
function is symmetric with respect to the left action of a certain subgroup of G, namely 
an isotropy subgroup of some element in the representation space of G. So, the next 
ingredient of our construction is a representation <&:GxVi— ►Vofa Lie group G in a 
linear space V; we denote it by 

&(g) -v for g G G, v G V. 

We denote also by <f> the corresponding representation of the Lie algebra g in V: 

d 



for £ G Q, v G V. 

e=0 



The map <j>* : g x V* i-» V* defined by 

is an anti-representation of the Lie algebra g in V * . We shall use also the bilinear operation 
o : V* x V h-> g* defined as follows: let u G V, y G V*, then 

(yov,$) = -{y,<f>(0 -v) V^Gg. 

(Notice that the pairings on the left-hand side and on the right-hand side of the latter 
equation are defined on different spaces. The operation o can be found in [19], where 
it plays the key role in the theory of the so called Clebsch representations; our notation 
follows II, 0|.) 

Fix an element p G V, and consider the isotropy subgroup G^ of p, i.e. 

G?M = {h : -p = p} C G. 

Suppose that the Lagrange function ~L(g,g) is invariant under the action of G^ on TG 
induced by left translations on G: 

L(hg,L h *g)=L(g,g), h€G®. 

The corresponding invariance property of LO(g,n) is expressed as: 

L®(hg,Sl)=L®(g,Sl), heG [p] . 

We want to reduce the Euler-Lagrange equations with respect to this left action. As a 
section {Gxg)/G^ we choose the set qxO p , where O p is the orbit of p under the action <£: 

O p = mg)-p, geG}cV. 

The reduction map is 

(g.fl)eGxg i-^ (tt,P) <E q x O p , where P = ^(g" 1 ) ■ p, 
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so that the reduced Lagrange function ■ X O p i— > R is defined as 
C {l) (n,P) = L {l) (g,n), where P = fcfo- 1 ) • p. 
The reduced Lagrangian £W(n,P) is well defined, because from 

P = c&(< 7l - 1 ) .p = d>( 5 - 1 ) . p 

there follows $ (c^si" 1 ) ■ P = P, so that e G [pl , and L®(#i,0) = L^(g 2 ,Q). 



Theorem 2.1. p], HI 



a) Under the left trivialization (g, g) \— * (g, O) and i/ie subsequent reduction (g, Q) i— > 
(f2,P), i/ie Euler- Lagrange equations ( \2. j\ ) become the following left Euler Poincare 
equations: 

f M = ad*0-M + V P £WoP, 

S ■ (2-4) 
1 P = -0(17) • P, 



where 



M = L* g n = VnC (l) . 



b) If the "Legendre transformation" 



is invertible, then (<L4) is a Hamiltonian system on q* x O p tritfi t/ie Hamilton function 

H(M, P) = (M, fl) - £ {l) (Q,P), 

with respect to the Poisson bracket given by 

{F 1 ,F 2 } = (M, [V M F 1 ,V M F 2 ]) 
+ (V P P 1 ,0(V A/ P 2 ) • P) - (V P P 2 ,0(V M Pi) • P) (2.5) 

for two arbitrary functions F\ j2 (M, P) : g* x O p *— > R. 



Remark 2.1. The formula ( |2.5| ) defines a Poisson bracket not only on g* x O p , but on all 
of g* x V. Rewriting this formula as 

{Fi,F 2 } = (M, [VmFx,V m F 2 ] } 
+ (P, f{V M F 2 ) • VpPi - <f (VmPi) • VpP 2 } 

one immediately identifies this bracket with the Lie-Poisson bracket of the semiproduct 
Lie algebra g x V* corresponding to the representation — <p* of g in V*. 

We shall also need a version of the above theorem for the case when the Lagrange 
function is invariant with respect to the right action of some isotropy subgroup G^ 
(A G V), i.e. 

L(gh,R h *g) = L(g,g), h€G^. 
In this case it is convenient to consider the right trivialization of the tangent bundle: 

(g,g) G TG i-> {g,u) G G x g, 
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where 

u = R g -i*g & g = R 9 *uj. 

Denote the Lagrange function pushed through the trivialization map by Jj^ r \g,uj) : G x 
g i-> R, so that 

LV(g,u)=L(g,g). 

The corresponding invariance property of the Lagrange function is expressed as 

LW(g/i,w) = LW( ?1 w), h€GW. 

The formulas for the reduction of the Euler-Lagrange equations are also slightly modified. 
As a section (G x g)/G^ we choose the set g x Oa, the reduction map being 

(g, u) G G x g i-> (u>, a) G g x Oa, where a = <&(g) ■ A , 

so that the reduced Lagrange function : g x Oa i— > R is defined as 

6 r \oj,a) = where a = • A 

There holds an analog of Theorem stating, in particular, that under the right trivial- 
ization (g,g) i— > (<7,u;) and the subsequent reduction (g,u) i— > (a;, a), the Euler-Lagrange 
equations ( |2.1|) become the following right Euler Poincare equations: 

fa = — ad* u ■ m — V a C^ o a, 

■ \ (2 - 6) 

a = <p(uj) ■ a, 

where 

m = R* g Tl = V U} &\ 



3 Lagrangian mechanics and Lagrangian reduction on GxG 

We now turn to the discrete time analog of these constructions. Our presentation of 
the general discrete time Lagrangian mechanics is an adaptation of the Moser-Veselov 



construction [37, 25| for the case when the basic manifold is a Lie group. The presentation 
of the discrete time Lagrangian reduction follows ||. Almost all constructions and 
results of the continuous time Lagrangian mechanics have their discrete time analogs. 
The only exception is the existence of the "energy" integral ( |2.2| ) . 

Let ~h{g,g) : G x G h- > R be a smooth function, called the (discrete time) Lagrange 
function. For an arbitrary sequence {g^ G G, k = ko, &o + 1> • • • > &i} ° ne can consider the 
action functional 

fci-l 

S=^L( 5fc , 9fe+1 ). (3.1) 
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Obviously, the sequences {g k } delivering extrema of this functional (in the class of varia- 
tions preserving g ko and g kl ) , satisfy with necessity the discrete Euler-Lagrange equations^. 

ViL(flfc,0 fc+ i) + V 2 L(g A ._i,5f fc ) = 0. (3.2) 

Here V{L(g, g) (V2^(g,g)) denotes the gradient of L(g, g) with respect to the first argu- 
ment g (resp. the second argument g). So, in our case, when G is a Lie group and not just 
a general smooth manifold, the equation (|3.2| ) is written in a coordinate free form, using 
the intrinsic notions of the Lie theory. As pointed out above, an invariant formulation 
of the Euler-Lagrange equations in the continuous time case is more sophisticated. This 
seems to underline the fundamental character of discrete Euler-Lagrange equations. 

The equation fl3.2|) is an implicit equation for g k +\. In general, it has more than 
one solution, and therefore defines a correspondence (multi- valued map) (g k -i,g k ) l— * 
(<7fc><7fc+i)- To discuss symplectic properties of this correspondence, one defines: 

U k = V 2 Ug k ^,g k )eT* k G. 

Then ( |3.2j ) may be rewritten as the following system: 

I II fc = -ViL(g k ,g k+1 ), ^ ^ 

I u k+1 = \7 2 Mgk,gk+i)- 

This system defines a (multivalued) map (g k ,H k ) i— > (g k+ i,H k+ i) of T*G into itself. More 
precisely, the first equation in ( |3.3[ ) is an implicit equation for g k +\, while the second 
one allows for the explicit and unique calculation of H k +ii knowing g k and g k +\ ■ As 



demonstrated in [37, |2f|, this map T*G h- » T*G is symplectic with respect to the standard 
symplectic structure on T*G. 

The tangent bundle TG does not appear in the discrete time context at all. On the 
contrary, the cotangent bundle T*G still plays an important role in the discrete time 
theory, as the phase space with the canonical invariant symplectic structure. The left 
trivialization of T*G is same as in the continuous time case: 

(g k ,U k )eT*G ^ (g k ,M k )eGx Q *, 



where 



Consider also the map 



M k = L* gk U k o U k = L* 1 M k . 



(fffc, ffifc+i) eGxGn (g k , W k ) £Gx G, (3.4) 



1 For the notations from the Lie groups theory used in this and subsequent sections see, e.g., M. In 
particular, for an arbitrary smooth function / : G ^ R, its right Lie derivative d'f and left Lie derivative 
df are functions from G into g* defined via the formulas 



{df(g),r,)= jj(^9) 



, (d'f(g),v)= ^/(ffe E ") 

e=0 ae 



and the gradient Vf(g) € T* G is defined as 

Vf(g) = R* g -rdf(g) = L* g -td'f(g). 
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where 

W k = 9^ 1 9k+i & Qk+i = 9kW k . 
Denote the Lagrange function pushed through ( |3.4[ ) by 

TL il) (g k ,W k ) = ~L(g k ,g k+1 ). 

Suppose that the Lagrange function L(<?, g) is invariant under the action of on 
G x G induced by left translations on G: 

Uhg,hg)=]L(g,g), heC®. 

The corresponding invariance property of L^(<7, W) is expressed as: 

L®(hg,W) =L®(jg,W), heG^. 

We want to reduce the Euler-Lagrange equations with respect to this left action. As a 
section (G x G)/G® we choose the set G x O p . The reduction map is 

{g, W) G G x G i ^ (W, P) € G x O p , where P = ^(g^ 1 ) ■ p, 

so that the reduced Lagrange function : G x O p i— ► R is defined as 

AW(W,P) =LW(c/,M/), where P = $( £ ,- 1 ).p. 

Theorem 3.1. [@, 1 

a) Under the left trivialization (g,g) *—> (g,W) and the subsequent reduction (g,W) h- > 
(W, P), i/ie Euler-Lagrange equations become the following left discrete time 

Euler Poincare equations: 



(3.5) 



f Ad* W^ 1 • M fc+1 = M fc + VpA^(W k , P k ) o P k , 
\ Pk+i = HW k l ) ■ P k , 

where 

M fc = 4AW(^- 1: P w )e fl *. 

b) If the "Legendre transformation" 

(W fc _i,i\_i) 6 G x O p m (M fc ,P fe ) G * x O pi 

_i, is invertible, then de/me a map (M fc , P fc ) h-> (M fe+1 , P fe+1 ) 
°/fl* x O p which is Poisson with respect to the Poisson bracket Ag. q) . 

We shall also need a version of the above theorem for the case when the Lagrange 
function is invariant with respect to the right action of some isotropy subgroup G^ 
(A e V), i.e. 

Ugh,gh)=Ug,g), hcG^. 

In this case it is convenient to consider the following map analogous to the right trivial- 
ization of the tangent bundle: 

(9k,9k+i) eG xG ^ (g k ,w k ) eGxG, 
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where 

w k = 9k+i9k l ^ 9k+i = Wk9k- 

Denote the Lagrange function pushed through this map by lj( r \gk, uu k ) ■ G x G *— > R, so 
that 

L^(g k ,w k ) = L,(g k ,g k+1 ). 
The corresponding invariance property of the Lagrange function is expressed as 

h^(gh,w) = U r \g,w), heG [A] . 

The formulas for the reduction of the Euler-Lagrange equations are also slightly modified. 
As a section (G x q)/G^ we choose the set q x Oa, the reduction map being 

(g,w) G G x G i-> (w,a) G G x Oa, where a = $(g) ■ A, 

so that the reduced Lagrange function A^ r ^ : G x Oa i— > R is defined as 

A^ r \w, a) = L,( r \g, w), where a = $(g) ■ A. 

There holds an analog of Theorem [O], stating, in particular, that under the right trivial- 
ization (g k ,gk+i) | — ► (9k, Wk) and the subsequent reduction (g k ,w k ) t— > («;&, a^), the Euler- 
Lagrange equations ( |3.3j ) become the following right discrete time Euler Poincare 
equations: 

Ad* w k ■ m k+ i =m k - V a A^(w k , a k ) o a k , 
afc+i = ®{w k ) ■ a k , 
where 

m k = d w A^ r \wk-i,a k -i) G 0*. 

The relation between the continuous time and the discrete time equations is established, 
if we set 

9k =9, 9k+i = 9 + £9 + O (e 2 ) , h(g k , g k+1 ) = eL{g, g) + O (e 2 ) ; 
P k = P, W k = I + eQ + 0(e 2 ), A^{W k , P k ) = e£.®(fl, P) + O (e 2 ) ; 
a fc = a, w fc = J + + O (e 2 ) , A (r) (w fc , a k ) = e£ (r) (u;, a) + O (e 2 ) . 

4 A multi-dimensional heavy top 

The basic Lie group relevant for our main examples is 

G = SO(n), so that g = so(n) 
(the "physical" rigid body corresponds to n = 3). The scalar product on q is defined as 

(f,»7) = -2 tr (^)' tV^S- 

This scalar product is used also to identify q* with q, so that the previous formula can be 
considered also as a pairing between the elements £ G g and i) £ g*. 
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The group G is a natural configuration space for problems related to the rotation of a 
rigid body. Indeed, if E(t) = (ei(i), . . . , e n (t)) stands for the time evolution of a certain 
orthonormal frame firmly attached to the rigid body (all e& G M n , (ejt, &j) = £jy), then 

E(t) =g- 1 (t)E(0) & e k (t)=g- 1 (t)e k (0) (1 < k < n), 

with some g(t) G G. The Lagrange function of an arbitrary rigid body rotating about a 
fixed point in a homogeneous field with a linear potential 

ip(x) = (p,x), 

is equal to 



H9,g) = -\tr(njn)-( P ,a), (4.1) 



where 



• Q = <7 1 g = g T g is the angular velocity of the rigid body in the body frame E(t); 

• J is a symmetric matrix (tensor of inertia of the rigid body); choosing the frame 
E(0) properly, we can assure this matrix to be diagonal, J = diag (Ji, . . . , J n ), which 
will be supposed from now on; 

• p G M. n is the constant gravity vector, calculated in the rest frame; 

• a = a{t) = g(t)A € R n is the vector pointing from the fixed point to the center of 
mass of the rigid body, calculated in the rest frame; A is the same vector calculated 
in the moving frame (where it is constant). 



Obviously, the function (4.1) may be rewritten also as 

L(g,g) = ±(n,J(n))-(P,A), (4.2) 
where J : so(n) i— > so(n) is the symmetric operator acting as 

j(n) = m + nj, 

and P(t) = g~ 1 (t)p is the gravity vector p calculated in the moving frame. The Lagrange 
function ([4.1[) is in the framework of Section 2, if the following identifications are made: 



• V = V* = M. n with the standard Euclidean scalar product; 

• The representation $ of G in V is defined as 

&(g) ■ v = gv for g G SO(n), »£] 

• Therefore the representation cp of q in V is given by 

<M£) ■ v = £v for ^ € so(n), v G M 
while the anti-representation (p* of q in V* is given by 

4>*{Q-y=-£,y for ^Gso(n), yG 
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Finally, the bilinear operation o : V* x V i— > g* is given by 
y ov = vy T — yv T = v A y for y,v £ 



The Lagrange function (|4.2| ) is manifestly invariant under the left action of the isotropy 
subgroup G^. Now Theorem [O] is applicable, which delivers the following equations of 
motion of a heavy top in the moving frame: 

f M=[M,S1] + AAP, 
\ P = -SIP, 

where 

M = V a C® = J (SI) O M jk = (Jj + j k )n jk . 



According to the general theory, the system (4.3) is Hamiltonian on the dual of the semidi- 



rect product Lie algebra e(n) = so(n) x R n , with the Hamilton function 

1 1 M 2 n 

H(M, P) = - (M, SI) + (P,A) = -J2 -r^r + E p ^ 



j<k 3 N k=l 

The correspondent invariant Poisson bracket reads: 

{Mi^Mkt} = M ie 5 jk - M kj 6 ei - M ik 5 je + M ej 6 ki , (4.4) 
{M ij ,P k } = P i 5 jk -P J 5 ik . (4.5) 

5 The multi-dimensional Lagrange top 

The multi- dimensional Lagrange top is characterized by the following data: J\ = J 2 = 
■ ■ ■ = J ra _i, which means that the body is rotationally symmetric with respect to the nth 
coordinate axis, and A\ = A2 = ■ ■ ■ = A n _\ = 0, which means that the fixed point lies on 
the symmetry axis. Choosing units properly, we may assume that 

Jl = J2 = ■■■ = Jn-l = f, J n = l-|, A= (0,0,...,0,1) T . 



The action of the operator J is given by 
M i:j = J{Sl)ij = | 
or in a more invariant fashion: 



aSlij, 1 < i,j < n — 1, 
Slij, i = n or j = n, 



M = J (SI) = aSl + (l-a) (SlAA T + AA T Sl) (5.1) 
= aSl-(l-a)AA(SlA). (5.2) 



Therefore, the formula ( [4.2j ) may be represented as 



L(g,g) = |<n,n) + L^(nA,SlA) - (P,A). (5.3) 
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The "Legendre transformation" (|5.1| ) 
there follows immediately 



is easily invertible. Notice that from ( |5.l| 



MA = UA, (5.4) 
and plugging this into the second term of the right-hand side of (|5.2| ), we find: 

CI = - M + — — A A (MA). 

a a 

The integrability of the multi-dimensional Lagrange top was demonstrated by Belya- 
ev P], who constructed "by hands" the necessary number of involutive integrals. A some- 
what easier proof is delivered by the Lax representation in the loop algebra sl(n+l)[A, A -1 ] 
twisted by the Cartan automorphism p0|, fjlf . 



Proposition 5.1. For the Lagrange top, the moving frame equations (4-5) are equivalent 
to the matrix equation 

L(X) = [L(X),U(X)], 



where 



L(A) 



M \A-\- l P 
XA T - A- X P T 



U(X) 



n xa 

XA T 



Proof. A direct calculation based on the formula (5.4). 



Integrability is not the only distinctive feature of the Lagrange top. Another one is the 
existence of a nice Euler-Poincare description not only in the moving frame, but also in 
the rest one. Rewriting (|5.3|) as 



^r(u,u) + ^—^-(u)a,u)a) - (p,a), 



(5.5) 



where u = gg' 1 , we observe that the Lagrange function of the Lagrange top is not only left- 
invariant (with respect to G^), but also right-invariant (with respect to G^ 4 '). Therefore, 
we may apply the formula (p.6[), which in the present setup reads: 



rh = [u, m] + V a £( r ) A a, 
a = ua, 



(5.6) 



Straightforward calculations based on ( |5.5| ) give: 

V a C {r) = -(l-a)to 2 a-p, 
m = Vw£ = olu + (1 — a) (ujaa7 + aa T u>) . 

The last formula implies, first, 

ma = acoa + (1 — a)uja(a, a) = ua, 

and, second, 

[u, m] — (1 — a) (oj 2 o) A a = [u, m — (1 — a) (ujaa T + aa T w)] = [u>, auj] = 0. 
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Plugging these results into (p.6|) , we finally arrive at the following nice system: 

{^ = aAP ' (5.7) 
I a = ma, 

where m = gMg~ l is the kinetic moment in the rest frame. This is a Hamiltonian system 
with respect to the (minus) Lie-Poisson bracket of e(n): 

{mij,m ke } = -mu8j k + m k jS ei + m ik dji - m ej 5 ki , (5.8) 

{rriij, cik} = -ai5jk + aj5i k , (5.9) 

with the Hamilton function 

H(m,a) = -(m,m) + (a,p). 



A remarkable feature of the system ( |5.7|) is its independence on the anisotropy parameter a. 

Proposition 5.2. The rest frame equations (\5.7Q of the Lagrange top are equivalent to 
the matrix equation 

i(X) = [£(\),u(\)], 

where 

m Xa — \~ 1 p 

Xa l - \- l p 




6 A discrete time analog of the Lagrange top 

An integrable discretization of the 3-dimensional Lagrange top was constructed in [0], 
where Lagrangian dynamics on G = SU(2) are used, leading to Euler-Poincare equations 
on su(2) x su(2). The multi-dimensional generalization presented here deals with La- 
grangians onG = SO(n), leading to Euler-Poincare equations on so(n) x M. n . This change 
of context (when compared with j?]]) required a modification of the kinetic energy terms. 
Of course, this is essential only in the discrete time setting, since in the continuous time 
situation the kinetic energy is given purely in terms of the Lie algebra q of the Lie group 
G, and the Lie algebras su(2) and so(3) are isomorphic. 

6.1 Rest frame formulation 

Consider the following discrete analog of the Lagrange function (|5.5|): 

H9k,9k+i) = -"trlog (21 + 10*. + w^ 1 ) 

- — — — log(l + (a k ,a k+1 )) - e{p,a k ). (6.1) 
e 

where w k , a k are defined as in Section 3: 

w k = fffc+ifffc 1 , a k = g k A, so that a k+ i = w k a k . 

The powers of e are introduced in a way assuring the correct asymptotics in the continuous 
limit e -> 0, namely h(g k ,g k+ i) w e~L(g,g), as g k = g and g k+ i w g + eg (see the end of 
Section 3). This is seen with the help of the following simple lemma. 
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Lemma 6.1. Let w(e) = I + euj + O (e 2 ) G SO(n) be a smooth curve, u G so(n). Then 

e 2 

tr log (21 + w(e) + w~ l (e)) = const - — (w, u) + O (e 3 ) . (6.2) 
For an arbitary a G W 1 : 

(a, w(e)a) = (a, a) — — (wo, uoa) + O (e 3 ) . (6.3) 
Proof. Let w = I + euj + e 2 v + O (e 3 ) . Then from ww T = I we get: 

u + v T + wu T = =4> v = -o; 2 + z;i, ui G so(n). (6-4) 

Hence 

2I + w + w T = 41 + eV + O (e 3 ) 
=> log (2/ + io + w T ) = const I + — lo 2 + O (e 3 ) , 
which yields ( |6.2| ). Similarly, we derive from ( |6.4[ ): 

/2 



T 

a toa = a 



"a + — a T w 2 a + O (e 3 ) , 



which implies (|6.3|). 



Writing in (|6.lD lOfcdfc for a^+i, we come to the Lagrange function A^ r \a k , w k ) = 
h(gk, gk+i) depending only on a k , w k , hence invariant with respect to the right action 
of G^ A K It may be reduced following the procedure of Section 3. 



Theorem 6.1. The Euler-Lagrange equations of motion for the Lagrange function (6.1) 
are equivalent to the following system: 

m k+ i = m k + ea k Ap, 

^ ^_ I+(e/2)m k+1 a (6-5) 
/ - (e/2)m k+1 



The map (m k ,a k ) i— > (m k+ i,a k+ i) is Poisson with respect to the bracket (\5.8j) , f |5. and 
the following function is its integral of motion: 

H e (m, a) = — (m, m) + (a,p) — -(m,p A a). (6-6) 



Proof. The right reduced Euler-Lagrange equations ( |3.6| ) have in the present setup the 
following form: 

J w^rrik+iWk = m k - a k A V a A (r \a k ,w k ), 
[ a k+ i = w k a k , 

where 

m k+ i = d w A^'(a k ,w k ). 
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To calculate the derivatives of A^ r \ we use the following formulas: 

d w trlog(2I + w k + w^) = -2^— -, (6.8) 
v K 1 Wk + I 

d w (a k ,w k a k ) = a k A (w k a k ), V a (a k ,w k a k ) = w k a k + wja k . (6.9) 
Indeed, the first one of these expressions follows from: 

(d w tr log (21 + w k + w - 1 ), r] ) = — tr log (21 + e en w k + 1 e~ eT >) 

d£ e=0 
= tr ((2/ + w k + w" 1 ) (w k - w^ 1 ) rj) = tr ((w k + I) _1 (i£; fc - I)rj) 

= -2((w k + I)- 1 (w k -I), V ). 



To prove the second one, proceed similarly: 

= (ak,V w ka k ) 



{d w (a k ,w k a k ),r)) = -^(a k ,e £V w k a k ) 
de 

tr (w k a k ajr]) = ^tr ((w k a k aj - a k ajwl) rj) = (a k A (w k a k ),r]) . 



6 = 



Finally, for the third expression we have: 

= (w k a k + wla k ,r]) . 



{V a (a k ,w k a k ),rj) = —(a k + er],w k (a k + er])) 
de 



e=0 



With the help of flO), (EM we find: 



2aw k -I 2(1 -a) a k A a k+1 

™k+i = — — r, (6.10) 

e w k + I e l + {a k ,a k+1 ) 



and 



w k 1 m k+1 w k + a k A V a A (r) (aAi, w k ) 

2a w k — I 2(1 — a) a k A afe + i 

= — z -T7 1 T ~ea k Ap = m k+ i - ea k A p. 

e w k + I e l + {a k ,a k+1 ) 

Comparing the latter formula with the first equation of motion in ( |6.7D , we find that it 
can be rewritten as 

m k = m k+1 - ea k Ap, 
which is the first equation in fl6.5|). 



To obtain the second one, derive from ( 6.10 ): 

e w k — I a k — a k+ i(a k ,a k+ i) 

-m k+ i a k+1 = a — a k+1 - (1 - a) — — — r , 

2 w k + I l + (a k ,a k+1 ) 

e w k -I a k {a k ,a k+ i) - a k+ i 

-m k+ i a k = a — a k - (1 - a) — —— r . 

2 w k + I ± + {a k ,a k+ i) 

Adding these two equations, we find: 

^-m k+1 (a k+ i +a k ) =a Wk \ (a k+ i + a k ) + (1 - a)(a k+1 - a k ) = a k+1 - a k . 
2 w k + 1 
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(On the last step we took into account that a k+ \ = w k a k .) This is nothing but the second 
equation of motion in fl6.5|) . 



The Poisson properties of the map ( |6.5| ) are assured by the version of Theorem 3.1 for 
right-invariant Lagrangians. 

It remains to demonstrate that the function fl6,6Q is indeed an integral of motion. This 
is done by the following derivation: 

H e (m k+ i,a k+1 ) = -(m k+ i,mk+i) + {a k +i - | m k +ia k+ i,p) 

= ^(m k +i,m k+ i) + (a k + ^m k+ ia k ,p) 

= ^(mk+i,m k+ i + epAa k ) + (a k ,p) 

= 7}( m k ~ epAa k ,m k ) + (a k ,p) = H e (m k ,a k ). 

(In this calculation we used the identity (ma,p) = (m,p A a) for m £ so(n), p,a £ M n ; 
notice that the scalar products on the both sides of this identity are defined on different 
spaces!) The theorem is proved. ■ 



Actually, the map (6.5) possesses not only the integral (^j) but a full set of involutive 
integrals necessary for the complete integrability. The most direct way to the proof of this 
statement is, as usual, through the Lax representation which lives, just as in the continuous 
time situation, in the loop algebra sl(n + 1)[A, A -1 ] twisted by the Cartan automorphism. 



Theorem 6.2. The map fl/5-3{ ) admits the following Lax representation: 

l k+ i(X) = v k 1 (\)£ k (\)v k (\), (6.11) 



with the matrices 



\ AW -A V / I — (e/2)ufc(A) 



m k Xb k -X 1 p\ ^ n _I+(6/2)a t (A) 

A- 

Xa k 

T 
k 



XaJ 



e \ e 2 



u k (X) = 

where the following abbreviation is used: 

b k =(l-~-m k )a k + j P . (6.12) 
Proof. Direct verification. Notice that it is most convenient to check ( |6.11 ) in the form 
(j + |u fc (A)) 4+i(A) (/ - |« fc (A)) = (/ - |« fc (A)) 4(A) (/ + |« fc (A)) , 

when it becomes polynomial in A. ■ 

This theorem provides us with a complete set of integrals of motion of our discrete time 
Lagrangian map: these are the coefficients of the characteristic polynomial det(4(A)— 
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Notice that these integrals of motion do not coincide with the integrals of motion of the 
continuous-time problem. To be more concrete, the integrals of motion of the map ( |6.5[ ) are 
obtained from the integrals of the continuous time problem by replacing a by b = a + 0(e) 
given in ( |6.12| ). As for the actual integration of our map in terms of theta-functions, we 



leave it as an open problem for the interested reader (cf., e.g., 28, [1 



6.2 Moving frame formulation 

It turns out that the equations of the discrete time Lagrange top in the moving frame for- 
mulation become a bit nicer under a little inessential change of the Lagrange function ( |6.lD , 
namely under replacing (p,a k ) in the last term on the right-hand side by (p,a k+ i). This 
modification does not influence the discrete action functional ( |3.1[ ), apart from the bound- 
ary terms. It is not difficult to see that this modification is equivalent to exchanging 
w k <-> u>£ 1 , a k <-> a&:+i7 which in turn is equivalent to considering the evolution backwards 
in time with the simultaneous change e i— > — e. Now express the discrete Lagrange function 
n|), with the above modification, in terms of P k = g k X p an d W k = 9 k 9k+i'- 

Hg k ,g k+ i) = ^ {l) (Pk, w k ) = ~ trlog {21 + w k + w^ 1 ) 

- 2(1 ~ a) log (1 + {A, W k A)) - e(P k , W k A). (6.13) 



Since W k = I + eSl + O (e 2 ) , we can apply Lemma 6J to see that 



A^(P k ,W k ) = eC^(P,n) + 0(e 2 ), 

where (P, Q) is the Lagrange function Q4.2| ) of the continuous time Lagrange top. Now, 
one can derive all results concerning the discrete time Lagrange top in the body frame 
from the ones in the rest frame by performing the change of frames so that 

M k = gl l m k g k , P k = g^p, A = g k x a k , 

and taking into account the modification mentioned above. Alternatively, one can derive 



them independently from and similarly to the rest frame results, applying Theorem 3.1 



the main result of which, the system (3J3), takes in the present setup the form 



f W k M k+1 W k l = M k + P k A V P A( l \W k , P k ), 
\ Pk+i = W^Pk, 

where 

M k+l = d{ /v A«\w k ,P k ). (6.15) 
Anyway, the corresponding results read: 

Theorem 6.3. The Euler-Lagrange equations for the Lagrange function ft6.13[ ) are equiv- 
alent to the following system: 

f M k+1 = W k l M k W k + eA A P k+1 , 
\ P k +l = W^Pk, 
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where the "angular velocity" W k € SO(n) is related to the "angular momentum" M k £ 
so(ra) via the "Legendre transformation": 



Ml. 



2a W k -I 2(1 - a) 
e 



A A (W k A) 



W k + I e I + (A, WkA}' (6,17) 

The map ( \6.1Q ), t \6.17 ) is Poisson with respect to the Poisson bracket ( 1^-4 A an d 

has a complete set of involutive integrals assuring its complete integrability. One of these 
integrals is 

H e (M,P) = ^(M,M) + (P,A) + ^(M,PAA). 

Proof. The derivation of ( |6.16| ), (6.17) is straightforward, like in Theorem |6.1| . Obviously, 
due to W k = I + eQ, + O (e 2 ), the equations of motion ( 6.16j ) approximate the continuous 
time ones ( |4.3[ ), while the "Legendre transformation" ( 6.17|) ap proximates (|5.2|). We dis- 
cuss now the inversion of the "Legendre transformation" ( |6.17 ). Obviously, it is trivially 
invertible if a = 1: then 



M, 



2 W k -I 



e W k + I 
For a general a we derive from ( 6.1 7| ) : 



I+(e/2)M k 
I-(e/2)M k 



w kA =L±sm^A. 



Indeed, we have: 



-M k A 
2 



a ■ 



W fc - J 
W k + I 

I 



e/2)M k 

A(A, W k A) - W k A 



(6.18) 



€ -M k W k A = a - 
2 W k + 1 



A -{I -a 



W k A - (1 - a) 



l + (A,W k A) ' 
A-W k A(A,W k A) 
1 + (A, W k A) 



Adding these two equations, we find: 



-M k (A + W k A) 



a 



W fc - J 
W k + I 



{A + W k A) - (1 - a){A - W k A) = W k A - A, 



which yields ( 6.18 ). Now plugging the expression for W k A through M k , A into the second 
term on the right-hand side of ( 6.17] ), we see that this equation is uniquely solvable for 
W k . Turning to the last statement of Theorem 5.3, we have: 



H e (M k+1 ,P k+1 ) = -(Mfc+i, Af fc+ i + eP k+1 AA) + (P k+1 ,A) 

= hw~ l M k W k - eP k+1 A A, W^M k W k ) + (P k+l ,A) 



= -(M k , M k ) + (P k+ i,A - - W~ l M k W k A) 

= l(M k , M k ) + (P k , W k A - e - M k W k A) 
\{M k ,M k ) + (P k ,A + e - M k A) = H e (M k ,P k ) 



(On the last but one step we used the formula Q6.18]) .) 
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We close this section with a Lax representation for the map ( |6.16D , ( 6.17] ). 



Theorem 6.4. The map ( 6. It ), ( 6.11 ) has the following Lax representation: 

L k+ i(X) = V k ~ 1 (X)L k (X)V k (X), 

with the matrices 

/ M k 



L k {\) 



Vk{X) ~ V 1 J I-(e/2)U(\y 
where the following abbreviation is used: 



XB k - X~ l P k 


U(X) = 



XA 
-XA T 



I+lM k ^A+ e -P k . 



7 The Clebsch case of the rigid body motion 
in an ideal fluid 



Clebsch [13] found an integrable case of the motion of a 3-dimensional rigid body in an 



ideal fluid, which was generalized in [26| for the n-dimensional situation. This problem is 
traditionally described by a Hamiltonian system on e*(n) with the Hamilton function 



1 1 n 



fe=l 



Here (M,P) G e*(n), so that M G so(n), P G W 1 , and C = {c jfc }™ fc=1 , B = diag(6 fc ) are 
some symmetric matrices. The equations of motion read: 



j M = [M,n\ + PA(BP), 
\ P = -£IP, 

where the matrix 17 G so(n) is defined by the formula 



(7.1) 



n 



CjkMj k . 



The "physical" Lagrangian formulation of this problem is dealing with a Lagrangian on 
the group E(n), left-invariant under the action of a whole group. However, it may be 
obtained also in the framework of Section 2, from a Lagrangian on SO(n), left-invariant 
under the action of an isotropy subgroup of some element p G M. n . These two different 
settings lead to formally identical results. 
So, one considers the Lagrange function 



L(g, g ) = c«\n,p) = \ (n,j(n)) + \ (p,bp). 



(7.2) 
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Here (g,g) € TSO(n), and £1 = g 1 g G so(n), P = g p 6 W 1 . It is supposed that 
J : so(n) i— > so(n) is a linear operator acting as 



The reduced equations of motion for this Lagrange function delivered by Theorem |2J. 
coincide with (7.1). 

The Clebsch case is characterized by the relations 

bi - bj bj -b k b k -bi 
+ — h 



C-jk C-ki 

which implies that 

_ bj ~ h 

Cjk — 

a,j — ak 



for some matrix A = diag (a^). The Lax representation found in [26] reads: 

L(A) = [L(A),C/(A)], 

where 

L(A) = XA + M + \- 1 PP T , U(X) = XB + n. 

In [34] we found an integrable Lagrangian discretization of the flow of the Clebsch 
system characterized by 

c Jk = T , T , h = J k , so that a k = J%. 
Jj + J k 

This flow may be considered as a particular case of the motion of a rigid body in a quadratic 
potential (notice that the kinetic energy term in Q7.2| ) is in this case typical for the heavy 
top, since J{Q) = JVL + Q J). Here we concentrate on another flow of the Clebsch system 
characterized by 

Cj k = 1, so that cifc = 
The corresponding Lagrange function reads: 

L(g, g) = (0, P) = \ <J2, fi> + \ (P, BP) , (7.3) 

so that 

M = fi, 



and the equations of motion become 



M = PA{BP), 
P = -MP. 



(7.4) 



The Hamilton function of this flow is given by 



H(M,P) = ±(M,M)-±(P,BP). 
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The Lax representation of this flow can be given in two equivalent forms: 

L(X) = [L(X),U + (X)] = [U^X),L(X)], 

where 

L(A) = XB + M + A _1 PP T , U+(X) = XB + M, C/_(A) = X- 1 PP T . 
Consider the following discrete time Lagrange function approximating (|7,3|): 

Ug k ,9k+i) = ^ l) (Pk,W k ) 
= ~ tr log (21 + W k + W^ 1 ) + * log + € -(P k , BP h fj . 

Here, as usual, W k = g^gk+i and P k = g~ x p. 

Theorem 7. 1 . The reduced Euler-Lagrange equations of motion for the Lagrange function 

Pk A (BP k ) 

RP„) ' 

(7.6) 



(7.5) 



([Oj read: 



M fe+1 = M fc + e 



l + (e 2 /4)(P fe ,PP fc 



J-(e/2)M fc+1 
n+1 -/+( e /2)M fe+i n - 

T/ie map (M k ,P k ) \— > (M k +i, P k +i) is Poisson with respect to the bracket (4-4)> 
This map is completely integrable and admits the following Lax representation: 



with the matrices 



L k {X) 



I+ ~, B 



L fc+ i(A) = V k (X)L k (X)V k 1 (X), 



where 



V k = (I + -M k )P k P£ (I--M k 



(7.7) 



l + (ey4)(P k ,BP k ) PkPk [I + - D 



Proof. Derivation of equations of motion is based on ( 6.14 ), ( |6.15| ). We have: 



+i 



2 Wife- J = J+(e/2)M fc+1 

e' W k + I k I- (e/2)Mfc+i 

which proves the second equation of motion in ( |7.6[ ), and 

ePPfc 



V P A^(W k ,P k ) 



l + (e 2 /4) (P k ,BP k y 
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which proves the first one. As for the Lax representation, it is verified by a direct calcu- 
lation. It is most convenient to check it in the form 

eA 



^-ft*) (AS + M k+1 + X^Pk+i) (i + ^Qfc) 
I + ^-^fc) ( XB + M k + A- 1 ^) (i - ^Qfc) 



(7.8) 



where 



1 ( . < 2 . • ., 



nk = TTWJWi^K) \ I+ -l B]PkPk 



Let us indicate the main steps of this calculation. Expanding (ffj]) in powers of A, we 
come to the following equations: 

A : M k+1 -M k = e(K k B-BQ k ); (7.9) 

A" 1 : Pfc+i - V k = - K k (M k+1 + M fc ) - | (M fc+1 + M fc )Q fc ; (7.10) 

A" 2 : (P fc+ i+P fc )Q fc -^(P fc+ i+Pfc) = |^ fc (M fc+1 -M fc )Q fc ; (7.11) 

A" 2 : ^ fc (P fc+ i-P fc )Qfc =0. (7.12) 

Here ( [7.9|) follows from the first equation in (|7.6| ). Equation ( [7.10 ) yields ( 7.12[) , which 
follows from the skew-symmetry of M k+ \, M k . To prove ( [7 . 1 0|) we notice that the first 
equation of motion in ([7.6[) easily yields the following expressions: 



£/, P,! ( / + \{M k+1 - A4)) , n k =(l- ~{M k+1 - A4)) P fe P fc T , (7.13) 



while the second equation of motion in (|7.6| ) shows that 

V k+1 = (/ - | M fc+ i) P fe P fe T (/ + | M fc+1 ) . (7.14) 

This is a companion formula to (f77| ). Now ( fTTOD follows directly from ( |7l3| ), 
(|7.7| ). Another consequence of these three formulas is: 

(p fe+ i + v k )Q k = n k {v k+1 + v k ) 

I - ~(M k+1 - M k j) P fc P fc T (i + -AM k+l - M k 



and this proves ( 7.11| ), since the right-hand side of this formula vanishes due to the skew- 
symmetry of M k+ i, M k . ■ 

We conclude this section with the following remarks. The discrete time system intro- 
duced here was found in the case n = 3 and in the su(2) framework by Adler M, who 
however did not use the Lagrangian formalism. It is well known that the restriction of 
the flow ( \l-4 ) to the symplectic leaf in e*(n) of the smallest possible dimension 2n is 
equivalent to the famous Neumann system which describes the motion of a particle on 
the surface of the sphere S = {i£ R n : {x, x) = 1} under the influence of the harmonic 
potential (1/2) (x, Bx). (The identification of the above mentioned symplectic leaf with 
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T*S is achieved via the formulas P = x, M = x Ap). It turns out that on this symplec- 
tic leaf our map ( |7.6| ) leads to Adler's discretization of the Neumann system [||, whose 
complete integrability was proved in [pq ]. The present construction delivers therefore a 
Lax representation for this discretization of the Neumann system - the problem left open 
in p§. 



8 Conclusion 

The models introduced in the present paper serve as further important examples of com- 
pletely integrable Lagrangian systems with a discrete time a la Moser-Veselov. The version 
of the discrete Lagrangian reduction leading to systems on duals to semidirect product 
Lie algebras was developed in Q. The list of relevant examples grows slowly, but now it 
already became quite representative: |2f| 0, and the present work. Several remarks 
are in order here. 

• Usually, discrete time Lagrangians lead to multi- valued maps (correspondences), 
cf. |25|]. This is also the case for the models introduced in |J7|, |34|] (discrete time 
Euler top and a discretization of the rigid body dynamics in a quadratic potential). 
Amasingly, the models of 0], as well as those of the present paper, lead to single- 
valued, and moreover explicit maps. Of course, one would like to be able to predict 
this outstanding property by just looking at the Lagrange function, but at present 
we do not know how to do this, since no deep reasons for this behavior are apparent. 

• A very intriguing and still not completely understood point is the capability of the 
discrete Lagrangian approach to produce completely integrable systems. It should 
be stressed that all the models listed above were found by guess, and the same holds 
for their Lax representations. There seems to exist (at least, at present) no regular 
procedure for finding decent integrable discretizations for finite-dimensional systems 
of the classical mechanics, like the rigid body dynamics. This is in a sharp contrast 
to the area of integrable differential- difference, or lattice systems, where the problem 
of integrable discretization may be solved in an almost algorithmic way, preserving 



the Lax matrix and the Hamiltonian properties, cf. [33]. 

For any completely integrable system, the involutive integrals of motion yield a set 
of commuting Hamiltonian flows (called a hierarchy in the infinite dimensional sit- 
uation). These flows share the Lax matrix with the original one, and there exist 
well-established procedures for finding the Lax representations for them. In the 
discrete time, the situation seems to be different. Although there is a recipe for 
producing commuting Poisson maps for a given one, provided its Lax representa- 
tion admits an r-matrix interpretation (which is the case for all the models listed 
above), there is no way, in general, for producing decent discretizations for concrete 
flows. Here "decent" means given by nice formulas in physical coordinates, or, if one 
prefers more precise terms, given by explicit Lagrange functions on Lie groups. In 
this connection it should be mentioned that in the only case when decent discretiza- 
tions are known for two different flows of the same hierarchy (the Clebsch problem, 
cf. and the present paper), these discretizations belong to different hierarchies 
and posess different Lax matrices and integrals of motion. 
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All in one, it should be said that while the Veselov's papers |37|, gg] appeared more than 
a decade ago and happened to give a strong impetus for the development of the whole 
subject of discrete integrable systems, the narrower area of integrable discrete Lagrangian 
systems pioneered there is still at the beginning of its development. 

Generally, we consider the discrete time Lagrangian mechanics as an important source 
of symplectic and, more general, Poisson maps. From some points of view the variational 
(Lagrangian) structure is even more fundamental and important than the Poisson (Hamil- 
tonian) one. (Cf. [18, 21], where a similar viewpoint is represented. Notice that the 
book [19 1 also dealing with the omnipresence of Hamiltonian systems on the semidirect 
product Lie algebras, takes an opposite viewpoint: the variational (Lagrangian) principles 
are derived there from the postulated Hamiltonian structures.) It would be important 
to continue the search for integrable Lagrangian discretizations of the known integrable 
systems. Also generalizations to the infinite dimensional case, e.g. to a discretization of 
ideal compressible fluids motion (see [18, [l~9[]), are highly desirable. 
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